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1 INTRODUCTION 

One of the consequence of Einstein General Relativity is that light 
is deflected by the gravitational field of an object with an angle of 



4GM 



where M is the mass of the deflecto r and h is the impact param- 
eter of the light ray teinstein igTTh . In the case of lensing of a 
background star by another star, this deflection creates extra im- 
ages of a background object where the images are too close to 
be resolved by the ground-base telescopes. This type of lensing is 
called gravitational microlensing. In 1936, Einstein said "there is a 
little chance of observing gravitational lensing caused by the stel- 
lar mass lenses". However, after several decades, using advanced 
cameras and computers, and following the propose of Paczynski 
(1986), the first microlensing event has been discovered by Alcock 
et al. (1993). 

The aim of the microlensing obs ervational groups 
such as MACH O. EROS and OGLE llXlcock et aL iOQOl ; 
lAfonso et al. 20031 ; IWvrzvkowski et al. 20091) has been to 
monitor stars the Large and the Small Magellanic Clouds in 
order to detect the massive compact halo objects (MACHOs) 
of the halo using gravitational microlensing. The passage of 
the MACHOs from the line of sight of the background stars 
causes an achromatic magnification of the background stars 
with a specific shape. Comparing the number of the observed 
events with that expected from the halo models have shown that 
MACHOs made up less th a n 20 per cen t of the Galac tic halo 
dMilsztain & Lasseire 200"TI ; iRahvar 20041 : iMoniez 2OI0l) . They 
can not cover all the dark matter budget of the halo. In this paper, 
we question further the existence of MACHOs on cosmological 
scales, and suggest a strategy for observation using the Hubble 
Space Telescope (HST). 



ABSTRACT 

In this work we propose a new strategy for the pixel-lensing observation of M87 in the 
Virgo cluster using the Hubble Space Telescope (HST). In contrast to the previous obser- 
vational strategy by Baltz et al., we show that a few days intensive observation with the 
duration of ~ 90 min in each HST orbit will increase the observational efficiency of the 
high-magnification events by more than one order of magnitude. We perform a Monte-Carlo 
simulation for this strategy and we show that the number of high magnification microlensing 
events will increase at the rate of 4.2 event per day with a typical transit time scale of 16 h. 
We also examine the possibility of observing mini-halo dark matter structures using pixel 
lensing. 

Key wor^^s.-gravitational lensing:Micro-darkmatter-galaxies:M87 

Gravitational microlensing towards the source stars beyond 
the Galactic halo can be performed with an observational method 
- so-called pixel lensing- which is different to the conventional 
method. When the background stars are too far and the field is too 
dense to be resolved, using this method a high magnification mi- 
crolensing event can change the light flux in each pixel. By record- 
ing the time variation of the lig ht in each p i xel, it is possible t o 
recognize a microlensing event l lCrotts 19921 : 1 Baillon et al. 1993h . 
An important factor for increasing the sensitivity of the observa- 
tions when using the pixel lensing is the size of point spread func- 
tion (PSF). Several groups as MEGA, AGAPE and ANGSTROM, 
have monitored M31. They have found few pixel-lensing candi- 
dates in this direct ion jCrotts & Tomanev I996l ; lAnsari et al. I999I : 
[kerins et al. 200^ . 

Gould (1995) proposed the extension of pixel lensing to cos- 
mological scales, using the observation by the HST towards the 
Virgo cluster. He proposed an investigation of the intercluster 
MACHOs and obtained a theoretical optical depth of r ~ 3 x 



(1) 



10 /, wher e / is the fraction of halo composed of MACHOs 
jGould I995I) . The detection threshold of an event was defined by 
the accumulation of the signal to noise of the events for the dura- 
tion that an event is observable. So, in this observation, one would 
expect to detect microlensing events with a long duration and a low 
peak in light curve. The estimation obtained for the rate of events 
with this observational threshold was F ~ 18/ day^^. Baltz et 
al (2004) used the HST for pixel lensing and monitored M87 for 
one month, using with the strategy of taking one date point per 
day. After data reduction, seven candidates remained. With further 
investigation of the light curves in two different colors they con- 
cluded that one candidate could be a microlensing event. Com- 
paring the observational result with the expectation from gener- 
ating synthetic light curves, they concluded that the fraction of 
MAC HOs in both the Mi lky Way and the Virgo halos are almost 
equal jAlcock et al. 2000l) . However, the problem with this statis- 
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tical analysis is tiiat one microlensing event has a large statistical 
uncertainty. In order to have a better estimation of the MACHOs in 
the Virgo cluster, we need to be observed more events. 

In this work we extend the work by Gould (1995) and Baltz 
et al. (2004), and we propose an alternative observational strategy 
using the HST. Our suggestion is to observe the very high magni- 
fication events in M87, which has a variation in its light curve in 
the order 1 d. In order to improve the sensitivity of observing of 
short duration events, we need to increase the sampling rate by the 
order of 1 h. We propose intensive short-duration observation of 
M87 of the order of 1 h using the Wide Field Planetary Camera 
3 (WFC3) of the HST. We suggest taking one date point in each 
orbit and collecting about 15 photometric data point per every 24 
h. This strategy will enable us to detect very high magnification 
short- transit pixel lensing events. 

First, we make a rough calculation to estimate the optical 
depth and the detection rate of the events using this observational 
strategy. In order to have a better estimation, we continue with a 
Monte-Carlo simulation for three observational programs of 1, 2 
and 3 d. We examine how sensitivity the results are to the duration 
of the observational program. We show that with 1, 2 and 3 d of in- 
tensive observation with the HST, it is possible to detect 4.2, 13.4 
and 18.5 events per day, respectively. We also simulate the obser- 
vation using the strategy of Baltz et al. (2004), with the cadence of 
1 d, and find 0.02 event per day for one month observations. 

The paper is organized as follows. In section[2l we introduce 
the Virgo cluster and obtain the theoretical optical depth and the 
event rate of the high magnification events from structures along 
the line of sight. In section[3] we perform a Monte-Carlo simulation 
to obtain the number of the high- magnification events found with 
the HST. The results are given in section|4] In section[5]we give our 
conclusion. 



2 PIXEL LENSING OF M87 

Virgo cluster contains ~ 1300 galaxies and is located about 16.5 
Mpc from us. It has am ap parent size of 8° and coordinate of RA~ 
12 hrs and DE~ 12° (Fouque et al. 200 ll ; iMei et al. 2007h . The 
brightest galaxy of this cluster is M87 which is located at the center 
of this structure. Similar to the Local Group, the mass of the Virgo 
cluster is mainly made up of dark matter. With microlensing of the 
Virgo cluster, it could be possible estimate the fraction of MACHOs 
in the halo of this structure. 

Because stars in the Virgo cluster cannot be resolved us- 
ing present ground- and space- based telescopes, the standard mi- 
crolensing technic cannot be used for the stars in the Virgo clus- 
ter. To make an estimation from the column density of the stars, 
we assume that galaxies in the Virgo cluster have the same sur- 
face density of stars as in the Milky way, (S — 50Mq /pc?). For 
HST with a PSF size of 0.067 arcsec in t he WFC3, th e number of 
stars inside the PSF would be about 1000 jPrive 2oTol) . To detect a 
microlensing event in a high blended background, we need a very 
high magnification event to change the flux of pixels inside the PSF. 
In what following, we provide a rough estimation from the optical 
depth and the event rate in terms of the transit time towards M87. 
In Section 3, we obtain a precise value of the event rate, based on a 
Monte Carlo simulation using /f 57 observations. 

The standard definition of the optical depth is the cumulative 
fraction of area covered by the Ein stein ring for all the distances 
from the observer to the source star jPaczviiski 19861) : 
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Figure 1. Number of stars inside the PSF size (solid line) indicated at the 
left hand side of y-axis and threshold for the impact parameter (dashed-line) 
indicated at the right hand side as a function of angular distance from the 
center of the galaxy. 
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where Re is Einstein radius, M is the typical mass of the lens, 
p{Di) is the mass density of the lens and Di is the distance between 
the observer and lens. 

Taking M87 as the target galaxy, the column density of stars at 
the center of this galaxy is much larger than that at the edge of the 
galaxy. Thus, the threshold for detecting a very high magnification 
event would depend on the position of the source star in M87. As- 
suming a spherical model for M87, we calculate the column density 
of the stars by integrating the number density of stars along the line 
of sight, S(r) = J n(r, z)dz. Then, the number density of stars 
for a given solid angle would be dN/dQ = D^T,. The number 
of stars inside the PSF is given by ANpsp ~ dN/dQ x Qpsp. 
Fig. m plots the number of stars inside the PSF in terms of sep- 
aration from the center of the galaxy in arcsec, using the galaxy 
model in McLaughlin (1999). The threshold impact parameter can 
be defined according to the blending inside the PSF. Because the 
threshold magnification should be more than the number of stars 
inside the PSF (i.e. Ath > ANpsp), the corresponding threshold 
for the impact parameter can be defined as Uth < (ANpsp)^^ . 
Fig. [T] shows the threshold impact parameter as a function of the 
angular separation from the center of the galaxy, as indicated by 
the dashed line. 

We can define a new optical depth according to the thresh- 
old impact parameter in such a way that Re = uthRp- Using the 
definition of the optical depth based on this reduced Einstein ra- 
dius, this relates to the conventional definition of optical depth as 
follows: 



f(ri) = Uth{n)T{n), 



(3) 



where, n is the direction of the observation and r is the standard 
definition of the optical in equation In what follows, we esti- 
mate the numerical value of the standard optical depth, resulting 
from the various structures along the line of sight. Finally, we ob- 
tain the optical depth of the observable high-magnification events 
by multiplying the optical depth by the threshold impact parameter, 
as in equation ([3). We note that the lenses along the line of sight can 
be in the halo of Milky Way, the disc of the Milky Way, the halo of 
the Virgo cluster, the halo of M87 or M87 itself. We calculate the 
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Figure 2. The optical depth towards the Virgo cluster in terms of the trans- 
verse angular separation from the center of the cluster for various structures. 
The solid line is the overall optical depth by accumulating the result from 
the different structures. 
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Figure 3. The reduced optical depth as a function of angular distance from 
the center of M87, resulting from the multiplication of the optical depth to 



optical depth towards all area of M87 and we provide the average 
optical depth of each structure, as follows. 

(i) We assume lenses in the galactic disc as the nearest struc- 
ture to the observer. The mass density for the lenses is given by 
teinnev & Tremaine 198% 



h 



(4) 



Here, h is tlie lengtli-scale of the disc, H is ttie tieiglit-scale and E 
is tlie column density of the disc at the position of the sun, Rq. 
Substituting p{R, z) into the equation (O and integrating over Di, 
the average value of the optical depth weighted by the number of 
background stars over M87, we obtain Ti ~ 1.4 x 10^'. 

(ii) For the Gal actic halo, we use NFW mass density profile 
jNavarro et al. 19961) 



ScPc 



(5) 



'{r/rs){l + r/rsy' 

Here, rs is the characteristic radius, pc — 3Ha /8itG is the present 
critical density, fh is the fraction of halo composed of MACHOs, 
r is the position of the lens from the center of galaxy and 5c = 
200c'' g{c) /3, where c is the concentration parameter and 

^^""^ " ln(l + c) - c/(l + c) ■ 

We set c = 18 and — lAkpc for galactic halo 

teattaglia et al. 200^ . We obtain the average optical depth from 
this structure m fh x 1.1 x 10~^, which is three order of mag- 
nitude larger than that of the Galactic disc. 

(iii) The third structure is the halo of Virgo cluster. We set 
the NFW parameters for t he Virgo halo as ~ 560fcpc and 
Scpc = 3.2 X lO^'^MM/pc^l lMcLaughlin 19991) . The optical depth 
is obtained as Tm ~ /u x 3.38 x 10~^, where /„ is the fraction of 
Virgo halo composed of MACHOs. We note that the optical depth 
from the Virgo cluster is also one order of magnitude larger than 
the Galactic halo. 

(iv) The lens can reside in the galactic halo of M87. We take 
the NFW model for the halo of this structure with the param- 
eters of rs = SOfc pe and c = 15 jGebhardt & Thomas 20091 ; 
iDohertv et al. 200^ . The optical depth resulting from this structure 
isr„ ~ /m87 X 7.62 x 10"^ 

(v) Finally, the lens can be in M87, so-called self-lensing 



events. We assume M87 is located at the center of the Virgo 
cluster with a spherical shape and a radius size of AOkpc 
dGebhardt & Tho mas 2009). For this structure, we use the mass 
density profile of (McLaughlin 1999): 



p{r) 



3 — 7 TbLb 



(7-4) 



(6) 



where 7 = 1.33, a = 5.1 ± 0.6 kpc. Lb = (5.55 ± 0.5) x 
10^" in, s,Ts = 14.6± 0.2 MqL, ' 



and r is the distance from 



the center of M87. The optical depth resulting from the self-lensing 
is-r„ ~ 1.73 X 10"^ 

We note that the halo of Virgo has the highest contribution 
in the optical depth towards M87. Fig. ^ shows the optical depth 
resulting from the various structures as a function of angular sep- 
aration from the center of the galaxy. Knowing the optical depth 
and the threshold impact parameter, we use the definition of the re- 
duced optical from equation ([3} and we obtain f in terms of angular 
distance from the center of M87 by multiplying the standard defini- 
tion of the optical depth in Fig.l|2]l by Uth{^) from Fig.[T] The result 
is shown in Fig[3] We summarize the result of the reduced optical 
depth in Table 1. Here, we assume that MACHOs contribute 20 per 
cent of the halo mass. 

We note that, in contrast to the standard definition of the op- 
tical depth, which is almost constant for all the directions towards 
the target structure, the reduced optical depth depends on the line 
of sight towards M87. It is small at the center of M87 and becomes 
larger at the edge of the structure. From the definition of the optical 
depth, we can obtain the rate for events for those with the impact 
factor smaller than Uth- From the dependence of the number of 
events on the optical depth, the rate of events is obtained as 

= ^-t, U) 

TobsNbg TV is 

Here, A^e is the number of the observed events, N^g is the number 
of the background stars during the observational time of Tabs and 
the reduced Einstein crossing time is defined as ts ~ uthts- Using 
the definition of r, the rate of events can be written in terms of the 
optical depth and the Einstein crossing time as 



2 r 

— —Uth- 

TV ts 



(8) 



TobsNbg 

One of the relevant time scales in the high-magnification microlens- 
ing events is which is defined as the full width at half 
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a b c d e overall 



T X 10" 


0.00074 


0.118 


3.622 


0.819 


0.0028 


4.563 




6.11 


5.39 


17.32 


9.02 


32.26 


15.85 



N^{l/day) 0.0032 0.607 5.794 2.513 0.013 8.93 



Table 1. Reduced optical depth (first row), transit time (second row) and the number of the events (third row) for each structure. Structures in each column are 
noted by (a) Galactic disk; (b) Galactic halo; (c) Virgo halo; (d) M87 halo; (e) self lensing and last column is mean value for the optical depth and transit time 
and the the overall number of detected events per day. 
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Figure 4. The event rate as a function of angular distance from the center of 
M87. The solid line is obtained from the rough calculation and the dashed 
line shows the result from the Monte-Carlo simulation. 
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Figure 5. The impact parameter normalized to the Einstein radius for 
micro-halos in terms of the impact parameter normalized to the size of 
micro-halo for the structures with various masses. 



maximum (FWHM) timescale of the high magnification events. For 
high magnification events, the magnification changes with the im- 
pact parameter as y4 = 1 /u. The transit time for these events relates 
to the Einstein crossing time as ~ 2y/3uotE- Hence, for each 
event, after fitting the observational light curve with the theoreti- 
cal light curve, we can obtain ti/2 in terms of ito and ts- Finally, 
taking the number density of background stars dN{6)bg / dO, as a 
function of angular distance from the center of the structure, we can 
obtain the differential rate of events as 

Tobsde dn tE{e) ' ^ ' 

Using this equation and the estimation for the time-scale of f^;, we 
obtain the rate of events as shown by the solid line in Fig.Q in 
terms of the angular distance from the center of M87. Integrating 
over 6 results in the number of the detectable events. Here, we ob- 
tain about ~ 9 high magnification microlensing events per day. We 
note that in this rough estimation, we did not use the quality of light 
curve on delectability of the events. Also we ignored the finite-size 
effect on the detections made with pixel lensing. In the section 3, 
we will re-analyse and simulate the light curves of M87 using the 
HST. Table (1) reports the reduced result of this preliminary anal- 
ysis and provide the reduced optical depth from the various struc- 
tures, ti/2 and and the rate of the events. 

Finally we exami ne the possibility of observing dark matter 
micro-halo structures dOoerdt et al. 20071 ; lErickcek & Law 20101) 
using pixel lensing. Taking the NFW profile for these structures, in 
order to have high- magnification events from these lenses, we have 
to decrease the impact parameter down to ~ 10~^. The problem 
with dark halo structures is that the mass of the structure enveloped 
by sphere with the radius of the impact parameter decreases to- 



wards the center of structure. Fig. ^ plots the impact parameter 
normalized to corresponding Einstein radius in term of the impact 
parameter normalized to the size of the structure. For microhalos, 
the normalized impact parameter, b/Rs is always larger than one 
and never reach below 10^'*, a favorable value for the pixel lens- 
ing. Our conclusion is that we can ignore possible contribution of 
these structures in pixel lensing. With current technology, it is not 
possible to detected the micro-halo structures. 



3 MONTE-CARLO SIMULATION 

In this section we perform a Monte-Carlo simulation to obtain the 
detection efficiency and the rate of microlensing events towards 
M87 using pixel lensing. In the first step, we simulate the stellar 
distribution in M87, using the color-magnitude profile of the stars 
of the galaxy. Then, we use the foreground structures to obtain the 
distribution of the lenses along the line of sight. In the pixel lensing, 
the magnification is suppressed by the blending and the finite size 
effect. We take into account these two effects in our simulation. 

To simulating the stellar distribution in M87, we use the 
theoretical Isochrones of the color-magnitude diagram diagram 
(CMD) obtained using the numerical simulation of the stel- 
lar evolution from the Padova model jMarigo et al. 200§) . The 
Isochrones are defined for the stars with different mass, color, 
magnitude but with the same matalicity and age. To simulat- 
ing the CMD, we select isochrones with the metalicity range of 
Z = 0.0004,0.008,0.004,0.001,0.019 and 0.030. For the evo- 
lution of the stars, we take the ages of the stars in the range of 
log{t/yr) = [6.6, 10.2] with interval of Alog(t) = 0.05 dex. 
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Figure 6. CM diagram generated by Padova Isochrones. The left panel is 
the absolute color-magnitude diagram and the light panel is the apparent 
CM generated according to the distance modulus, extinction and reddening 
of the stars. We show a sample of microlensed stars with a high signal to 
noise of Q = 30 (cross sign) and low signal to noise of (plus sign) (3 = 5 
in the Pixel lensing. 



For each age range, we take a metalicity according to the age- 
metalicity relation, and the distribution function for the metahcity 
from Twarog (1980). In order to generate the CMD of the stars at 
the present time, we need an initial mass function of the stars and 
the star formation rate function for M 87. The initial m ass function 
is chosen from the Salpeter function dSalpeter 1991) . For the star 
formation rate, we assumed a positive rate dN/dt > similar to 
that in our Galaxy dCignoni et al. 20061 : iKarami lOld) . The result 
of our simulation in the absolute CMD is shown in the left panel of 
Fig.©. 

In order to generate the apparent color-magnitude of the stars, 
we calculate the distance modulus of each star in M87 and add 
it to the corresponding extinction. Because dust in the interstellar 
medium is proportional to the Hydrogen density, we can use the 
column density of hydrogen {Nh) along the line of sight as an in- 
dicator of amount of exti nction. For the V and I bands , the amount 
of extinction is given by jWeingartner & Draine 200ll) : 



Av{mag) = 5.3 x 10 NH{atom/cm^), 
Ai{mag) = 2.6 x 10'^'^ Nh {atom/ cm'^). 



(10) 



For the distribution of the hydrogen in M87, we choose the theoret- 
ical distribution fu nction of hy drogen, which has the best fit to the 
X-ray observation ( iTsai 19941) : 



ne{r) = no 



(r/ai)-° 
1 + (r/ai)' 



(11) 



' and 



Here, ai = 9.51fcpc, a = 0.36, no = 3.51 x 10~^ atom/cm" 
r is the distance from the center of the structure. 

Finally, the apparent magnitude and color of each star are 
given by 



my 



Cob 



Mv + 5 log 



d{pc) 



10 



+ Av, 



Co + Ai 



A,, 



(12) 



Here, Ca is the absolute color. My is the absolute magnitude in 
the visible band and d is distance of stars in parsec. We plot the 
apparent color-magnitude of the stars of M87 in the right panel of 



Fig.([6ll, including the effect of the distance modulus, extinction and 
the reddening. 

In the next step, we generate the light curve of the microlens- 
ing events by choosing the parameters of the lens as the mass, the 
velocity and the location, according to the corresponding distribu- 
tion functions. The location of the lenses from the observer is calcu- 
lated according to the probability function for the location of lenses 



dV/dx oc p{x) \/ x(l — x) 

, where x = Di/Ds and p{x) is the mass density profile of 
the lenses along the line of sight. The mass of the lenses, in 
units of solar mass , are chosen f rom the Kroupa mass function 
jKroupa et al. 19931 : iKroupa 200 ih . f oc M~" , where a = 0.3 
for 0.01 4: Ml < 0.08, a = 1.3 for 0.08 ^ Mi < 0.5 and 
a — 2.35 for 0.5 ^ Mi ^ 1.0. The velocity of lenses in the halo 
is chosen from the Maxwell-Boltzmann distribution with the dis- 
persion velocity indicated by the Virial theorem. For the Galactic 
halo, the Virgo halo and the M87 halo, the dispersion velocities are 
set as a ~ 156 fcms~^, 1000 kms^^ and 350 km s^^, respec- 
tively jPehnen et al. 20061 : lOebhardt & Thomas 20091) . The distri- 
bution for velocity of stars in the Galactic disc is also chosen from 
ellipsoid of the dispersion velocity jBinnev & Tremaine 1987h . The 
velocity of stars in M87 is chosen from Maxwell-Boltzmann func- 
tion with a ~ 360fc?ns^^ jRomanowskv & Kochanek" 200 ih . The 
minimum impact parameter uo in this simulation has to be chosen 
geometrically uniform in the range of uq G [0, 1]. However, in the 
Monte-Carlo simulation, we generate a uniform impact factor in the 
logarithmic scale to have enough statistics of the high- magnifica- 
tion events and to reduce the computation time. In order to recover 
the uniform distribution of impact parameter, we consider a weight 
proportional to the impact parameter uq when counting the number 
of events. 

To calculate the finite-size effect, we need the size of the stars 
projected on the lens plane. For the main-sequence source stars, the 
radius of the stars is obtained according to the mass-radius relation 



(13) 



Here, the parameters are normalized to the Sun's values. Red giants 
(RGs) also follow the relation of dHavashi et al. 1962i) 



^RG^^RG ~ const. 



(14) 



The formalism of the fin ite-siz e effect in the gravitational lensing 
was obtained by M aeder |l973h and developed later for microlens- 
ing by Witt & Mao ( Il994l) . The relevant parameter in the finite size 
effect is p^/uo, where is the projected radius of the star on the 
lens plane normalized to the Einstein radius. For the population of 
the source stars in M87 and lenses in various structures along the 
line of sight. Table l|2j gives the average value of p*/ito, with the 
corresponding p^, and uq for the observed events in each structure. 
We examine the finite-size effect on the main sequence and RG 
stars. We find that 43% per cent of stars with an absolute magni- 
tude less than zero are RG, but for the microlensed stars in this 
range only 17 per cent are RGs. 

In the pixel lensing, because of the high blending in each pixel, 
the signal to noise ratio from microlensing must be so high to be 
distinguishable from the background brightness of the host galaxy. 
For a background star with the apparent magnitude of m, the extra 
photons during the magnification of the star are obtained by 

-0.4{m-m.'P) 



SN ^ {A - 1) X t^^p X 10 



(15) 



Here, A is magnification factor of the source star, texp is exposure 
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a b c d e overall 
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T X 10'' 


0.00139 


0.219 


6.758 


1.525 


1.734 


10.238 


2 


Ne{l/day) 


1.2e-6 


0.0022 


4.2 


0.018 


0.0072 


4.2 






a 


b 


c 
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e 


average 


3 


P*/«0 


0.002 
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Table 2. Result from the simulation of one day observation by the HST with the WFC3 camera. The columns are (a) Galactic disk;(b) Galactic halo; (c) Virgo 
halo; (d) M87 halo; (e) self lensing. The first and second row represents the optical depth and the number of the events resulting from each structure and the 
last column indicates the overall numbers. The third line is the mean value of /mq for the observed events in each structure. The forth row represents the 
average transit time scale of the events, the fifth row is the average of maximum magnification, the sixth row is the impact parameter and the seventh row is 
the apparent magnitude of the source stars. The last column represents the average value for the parameters of the lenses. 



time and m^^ is the zero-point magnitude equivalent to the flux 
= 1 photoelectron"^, which is equal to 25.10 and 26.08 mag in 
/-band (F814W) and V-band (F606W), respectively, for the WFC3 
camera of the HST. We note that the observation by Baltz et al. 
(2004) were performed with the WFC2, which had a lower zero- 
point magnitude of 23.86 mag for the / band. The exposure time 
can be obtained from the visibility of M87 by HST. For M87, which 
is located at DEC = 12°, the visibility would be about 52 min in 
each HST orbit. We use the /-band (F814W) filter for observation 
in the simulation. We show that the background events, such as 
novas, cannot pass our light curve criterion for a short observation 
period. Thus, we do not use another filter, such as the R band, for 
the variable identifer and devote all the observation time in the / 
band to increase the signal-to-noise ratio. Within each orbit, we 
assume six 260 s exposure with an overall exposure of 1560 s in 
/ band. The noise is considered as the intrinsic Poisson fluctuation 
inside the PSF. This is proportional to the square of the number of 
the photons received from the background sky and the host galaxy: 



(16) 



Here, Q.psF i& the characteristic area of the PSF, ^sky is the bright- 
ness of the background sky (set to 21.(Smag as^'^ and /i is the sur- 
face brightness profile of M87. In our Monte-Carlo simulation, we 
accept the events with the following criterion: 



Q = /^(5iV>/\/iV)2 > 



(17) 



Here, the s ummation is perf ormed over the observed points of the 
light curve dOould A. 19961) and Q^r it is the critical value for the 
signal-to-noise of the detection (chosen to be Qr ~ 5 for a relaxed 
criterion and Qs — 30 for a strict criterion). The is equivalent 
to the Xbi fit of the light curve to the baseline. 

As the second criterion, we adopt those events in which the 
peak of the light curve resides within the duration of the HST ob- 
servation, (i.e. to £ [O.Tobs]^ means that we detect events with 
both rising and declining parts of the light curve). For the short 
term events, the probability of fulfilling this criterion is higher 
than for the long duration events. We add another criterion, which 
is sensitive to the fast variation of the light curve, for which we 
detect at least three consecutive data points, 3a above the back- 
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Figure 7. The theoretical light curve of a high magnification microlensing 
event including the blending of the background stars with the simulation of 
the observation with the cadence of ~ 90 min (top panel data points). In 
the bottom panel we show the same event for the duration of 40 days. 



ground flux. This means that we have, at least, a 3 x 3(t devia- 
tion of the microlensing light curve from the bassline fit. In other 
words, the difference between the two different fits should be at 
least Ax' = xli " xL > 27. 

In Fig.l|6ll, we identify the observed events that result from the 
simulation in the CMD using '+' symbols for the relaxed criterion 
and ' X ' symbols for the strict criterion. We also plot a typical sim- 
ulated light curve for a high-magnification microlensing event in 
Fig|7] The error bars are la results from the Poisson fluctuations in 
the number of photons received by the detector, according to equa- 
tion il6l . shifted by a Gaussian function from the theoretical light 
curve. The top panel of the Fig. |7] represents the peak of the theo- 
retical light curve and the simulated data points with the duration 
of ~ 90 min. The whole light curve is shown in the bottom panel. 
Similar to the distribution ofts towards a given source, we plot the 
distribution of the relevant transit time of ti/2 in Fig. [8] The overall 
number of the events are normalized according to Table [2] 

One of the backgrounds that we should take into account in 
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Figure 8. The distribution function of the observed events witli tlie 

criterion of Q > 5. 

the simulation of the light curves is the hot pixels of the CCD re- 
sulting from the cosmic rays. The probabilit y of a hot pixel in th e 
HST foi 1000 s exposure time is 1.5 ~ 3% jSirianni et al. 2005h . 
We take into account the effect of cosmic rays and discard the hot 
pixels from the light curve. The other backgrounds that can be mis- 
taken for a microlensing events are the variable stars and novas. Th e 
annual average Novas in M87 is about ~ 64 ( [Madrid et al. 20071) . 
So, during the observations over 1 or 2 d, we can expected one 
event, at most, to be detected. Also, the variable stars do not have 
enough variation in the light curve to be detected in the high blend- 
ing medium. We ignore the effects of novas and variable stars. 



4 RESULT 

In this section we report the results of the Monte-Carlo simulation 
as the number of events and the mean value for the parameters of 
the lenses. For each structure along the line of sight, we define a 
threshold form the simulation where a microlensing event is de- 
tected in the range of < Uth- The overall number of events can 
be obtained by summing the number of events for each structure as 
follows: 

Here, Ns is the number of the structures, the superscript 'i' cor- 
responds to the ith structure along the line of sight and e'its) 
corresponds to the detection efficiency. < t'^/e^{tE) > is the av- 
erage value for the ratio of the reduced Einstein crossing time to 
the detection efficiency with the threshold impact parameter of uth 
for each structure. The integration is performed over the surface of 
M87 using the column density of the stars from the model. We also 
obtain the differential number of events dN /dO as a function of 
angular separation from the center of the structure, as represented 



Figure 9. The detection efficiency as function of the log{uo) (top panel) 
and log(ti/2/'') (low panel) for two strategies of short cadence and long 
cadence observations. The short cadence observation with one day, two days 
and tluee days ai'e shown with green dotted line, blue dots-dashed line and 
black solid line, respectively. The efficiency for the long cadence observa- 
tion with the duration of one month is shown with red dashed line. 

by the dashed line in Fig|4] Unlike the results from the rough esti- 
mation, the results from the Monte Carlo simulation show that the 
central part of the structure is better for pixel-lensing. Our simu- 
lation results in 4.2 pixel lensing for 1 d of observations with the 
HST. The corresponding parameters of the events are shown in Ta- 
ble |2] In order to test how sensitive the detection is to the duration 
of the observation, we run the Monte-Carlo simulation for the du- 
ration of 1, 2 and 3 d using both with the new strategy and that of 
Baltz et al. (2004) for one month. We report the daily numbers of 
events in Table l|3). 

The results for the observational efficiency is shown in Fig 
|9l where we plot the detection efficiency in terms of log{uo) (top 
panel) and log{ti/2) (bottom panel) for two observational strate- 
gies and four observational durations. The minimum impact pa- 
rameter is limited by the finite-size effect below log{uo) — —2.5, 
where the efficiency remains constant. The efficiency diagrams pre- 
dict that intensive observations using the HST over one or a few 
days, with the duration of the orbital time, is much more efficient 
than observations over one month, taking data point per day. The 
new strategy increases the observational efficiency by almost one 
order of magnitude. Table |3] summarize the results of simulation 
for various observational strategies. 



5 CONCLUSION 

Gould (1995) proposed an observational strategy for pixel lensing 
of the Virgo cluster using the HST, taking one data point per day 
for the duration of one month. Baltz et al. (2004) performed the 
observations with the WFC2 camera of the HST and found seven 
variable candidates, from which one was shown to be a microlens- 
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Table 3. The first three columns are for one, two and three days observation. The forth column is the result of simulation for the observation with the duration 
of 30 days with the strategy of taldng 30 data points during one month. All the simulation is done based on using WFC3 HST camera, except the last one 
with WFC2. The first row shows the number of events per one day. The second row indicates the average transit time of the events in hour and the third row 
indicates the average Einstein crossing time in day. 



ing candidate. They used one microlensing candidate and put limit 
on the contribution of MACHOs in the Virgo halo. 

In this paper, we have proposed a new observational strategy, 
emphasizing the detection of the shor-duration very high mag- 
nification events. To detect such events, we need high cadence 
observation by HST, which means one observation per orbit for 
the duration of one to few days. First, we made a rough estimation 
of the number of observed events using this strategy. For the 
detailed study, we performed a Monte-Carlo simulation to obtain 
the number of pixel lensings events. With HST observation over 1 
d, and considering 20 per cent of the Virgo halo to be composed of 
MACHOs, we expect to see ~ 4.2 events. If we look at the result 
of the observations over 2 or 3 d, there are more than two and 
three times the number of the events than for observation taken 
over 1 d. We note that, with the new strategy, taking observations 
over 15 orbits in 1 d uses half of the observational time than one 
month of observations with the former observational strategy, and 
it increases the number of events substantially. This new strategy 
could provide new possibility for pixel-lensing observations on 
cosmological scales. 
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